Crucial role of decoherence for electronic transport in molecular wires: Polyaniline as 

a case study. 
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In this work we attempt to elucidate the nature of conductivity in polymers by taking the acid- 
base doped polyaniline (PAni) polymer. We evaluate the PAni conductance by using realistic ab 
initio parameters and including decoherent processes within the minimal parametrization model 
of D'Amato-Pastawski. In contrast to general wisdom, which associates the conducting state with 
coherent propagation in a periodic polaronic lattice, we show that decoherence can account for 
high conductance in the strongly disordered bipolaronic lattice. Hence, according to our results, 
there is no need of considering a mix model of "conducting" polaronic lattice islands separated 
by "insulating" bipolaronic lattice strands as is usually assumed for PAni. We find that without 
dephasing events, even very short strands of bipolaronic lattices are not able to sustain electronic 
transport. We also include a discussion of specific mechanisms that should be involved in decoherence 
rates of PAni and relate them with Marcus-Hush theory of electron transfer. 

PACS numbers: 72.80.Le, 72.10.-d, 73.20.Jc, 72.15.Rn 



I. INTRODUCTION 

In the late 70's, Alan MacDiarmid, Alan Heeger and 
Hideki Shirakawa led the investigations which put con- 
ducting polymers at the center stage by unraveling the 
transition from insulator to metal upon doping of poly- 
acetilene The following decades, these materials 

encountered numerous technological applications 0, Q. 
The novelty of polyacetilene's physical properties, e.g. 
transport through solitonic excitations Q, made it the 
most intensively studied conductive polymer. However, 
the interest then shifted to polyanilines and related 
compounds because they are inexpensive, stable and easy 
to made. 

In spite of its long history, polyaniline (PAni) became a 
new paradigm for polymeric conductors as it shows a dra- 
matic increase in conductivity either by acidic treatment 
or by electrochemical oxidation. In spite of this fact, 
the physical basis of its transport mechanism and of the 
insulator-metal transition proved more elusive. Starting 
from a semiconducting PAni in an emeraldine base form 
(Fig. HJ-a), protonation leads to an internal redox reac- 
tion that converts it into a metal (emeraldine salt). In 
order to account for the highly conducting nature of this 
doped polymer there are two well established models that 
imply two different lattice arrangements. These are asso- 
ciated to the appearance of two possible charged defects 
upon protonation. On one hand, the polaronic lattice 
(PL), which describes a lattice of Nitrogen bridged ben- 
zene rings that becomes fully periodic in the case of 100% 
of protonation. Even when one of every two Nitrogens 
is in the form supporting a polaron, the correspond- 
ing -pz electrons form a collective band of Bloch extended 
states which, being half-filled, behaves as a metal (Fig. [T]- 
c). On the other hand, in a crystalline bipolaronic lattice 
(BL) the protonated quinoid units {NH+ = Q = NH+) 



are bridged by three benzene rings. The electron tun- 
neling between neighbor NH^ units, leads to a bonding 
basic unit that justifies a bipolaronic description. Natu- 
ral disorder appears through the fluctuation of the bridge 
length (Fig. [T]-b). Hence, while further tunneling be- 
tween units could be possible, within the standard wis- 
dom, disorder ensures localized eigenstates that prevent 
propagation^. Galvao et al. [?[ concluded that disor- 
der should be an essential ingredient in these systems. 
They made molecular orbitals calculations of the elec- 
tronic structure of PAni chains which showed that dis- 
order pulls the Fermi energy down through the localized 
states of the valence band. Later on, Wu and Phillips 
l^i agreed with Galvao in the role of the protonation, 
further showing that induced disorder can be identified 
with a Random Dimer Model (RDM) By adopt- 

Landauer's view that "conductance is transmission" 
the current motto of molecular electronics 12], it 
was proved that the short range order of the RDM pro- 
duces a set of delocalized or propagating states [9l. [l0lfl3| . 
This opened the possibility that Fermi energy might lie 
in a delocalized region. However, Farchioni et al. [14|. by 
using an ab initio parametrization, made a detailed tight- 
binding based study of PAni-HCl comparing the BL and 
PL models [l5|. They showed that even when the BL 
model exhibits extended states, its Fermi energy is far 
from the high transmission regions. These ideas seemed 
to support the PL as the only PAni emeraldine salt capa- 
ble of metallic behavior. Indeed, the observation of Pauli 
susceptibility on conducting samples was attributed to 
extended states in a polaronic lattice [l6|. However, first- 
principle energy stability calculations point into the op- 
posite direction. A BL is by far the more stable energy 
configuration when compared to a PL [Tt} or its vari- 
ants [l^ . A picture that could unify these conclusions is 
that of segregated metallic (PL) regions and insulating 



(BL) domains. Transport would be mediated by hop- 
ping between metallic fibers in the polymer backbone 
|19l |. However, it is not clear that such structures could 
give a lower free energy than a pure BL one. Besides, in 
this model disordered islands would constitute the con- 
ductance bottle-neck for which a microscopic description 
is lacking. Further emphasizing the role of BL, it was 
recently suggested that susceptibility experiments could 
not be used to rule out the bipolaronic structure from 
conducting samples. This is because an internal chemi- 
cal redox equilibrium between bipolaronic structures and 
a number of polaronic defects with Curie susceptibility, 
should manifest as an overall susceptibility whose tem- 
perature dependence would be indistinguishable from the 
Pauli paramagnetism [20*1. In summary, the early works 
associate the conducting state of PAni with periodic or- 
der because the existence of extended Bloch eigenstates 
is a condition for coherent propagation. 

In order to account for the surprising frequency depen- 
dence of the dielectric constant and of the conductivity 
observed on conjugated polymers [2l| . Prigodin and Ep- 
stein i22j suggested a new mechanism of charge trans- 
port. They argued that the metallic state of polymers 
like PAni is sustained by a granular picture of transport 
where metallic islands, separated by amorphous material, 
interact through intrachain resonant tunneling events in 
a quasi 1-D variable range hopping theory. However, 
after an energy scale analysis. Martens et al. [1^ H^l 
arrived to the conclusion that intrachain charge carrier 
delocalization should extend over several grains. In con- 
sequence, there is some critical mechanism that governs 
the formation of truly delocalized states. They pro- 
pose a quasi 1-D model of weakly coupled disordered 
chains with phase-breaking events that are modeled in 
the Landauer-Biittiker framework. In this case the 1-D 
Schrodinger wave function picture for a single chain re- 
mains essentially correct with the additions of a finite 
lifetime, i.e. decoherence, due to dephasing events. Their 
source can be multiple, ranging from electron-phonon 
coupling [H,!!^ to through-space electron-electron inter- 
actions between charge fluctuations (27j , or even a weak 
interchain coupling [28|. Increasing the interchain cou- 
pling eventually will give rise to a transition from a quasi 
1-D to a fully 3-D behavior as demonstrated by numerical 
simulations [H, Is^]. However, for conjugated polymers 
such as PAni, the interchain charge transfer is weak and 
a 1-D model that includes decoherence should be a good 
approximation. Within this framework. Martens et al. 
invoke dimensional arguments that explain the anoma- 
lous frequency dependence of the dielectric constant and 
conductivity of several polymers. However, their con- 
clusions are based on estimations of the relevant system 
quantities. 

In this work we attempt to elucidate in detail the na- 
ture of conductivity of polymers by taking the PAni bipo- 
laron lattice structure as the case study. For that pur- 
pose, we use realistic ah initio based tight-binding pa- 
rameters which can be easily reduced jsii] to the min- 
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FIG. 1. Three different forms of PAni: (a) Emeraldine base, 
and two lattice models after a doping process, (b) Bipolaron, 
(c) Polaron 



imal parametrization of the D'Amato-Pastawski (DP) 
model [32] . This provides a simple solution to the other- 
wise complex Keldysh formulation of transport [13, 113] . 
Indeed, this strategy was applied before to PAni by 
Maschke et al. and Schreiber et al. (33 - [36j . How- 
ever, they focused on already conducting PAni polaronic 
chains that are affected by decoherence and/or interchain 
coupling. This also sums to the recent efforts in includ- 
ing decoherent processes in molecular electronics [37l - [39| . 
At this point a comment on the meaning of decoherence 
seems in order since nowadays this term is mostly used 
in connection with quantum information theory, where 
it denotes the (full or partial) collapse of a pure quan- 
tum state due to usually unspecified interactions with 
the environment. This might seem different from the 
usual language of electronic transport in the solid state. 
There, one is used to deal with specific interactions, such 
as those with phonons, magnetic impurities or other elec- 
trons, which may involve transitions with given selection 
rules, e.g. |fc) — > \k + q). These transition probabili- 
ties are evaluated within a Fermi Golden Rule (FGR). 
In this approximation, the coupling with the extra de- 
grees of freedom, the environment, prevents the interfer- 
ence among the component remaining in \k) with that 
in \k -\- q). Thus, with some probability the environment 
"measures" an electron in the state \k) and "re-injects" 
it incoherently in \k + q). These processes, which usually 
involve some inelasticity, are inherently different from 
the elastic scattering with imperfections and impurities 
that produce the interferences leading to localization @ . 
Quite often one realizes that, regardless of the specific se- 
lection rule, the relevant role of interactions is just to pro- 
vide for decoherence, a mechanism that competes with 
the coherent scattering that results in localization. This 
is precisely the spirit of the "local phonon" bath or the 
fictitious voltage probes that lead to the imaginary site 
energies introduced in the DP model [H, H^]. Thus, the 
idea is that decoherence from the system-environment in- 
teraction provides a knob that sweeps transport between 
a Mott's variable range hopping regime and a typical 1- 
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D metal [40l . |4l| . Transport in conducting PAni should 
occur somewhere at this crossover. 

We start in section II by recalling the Landauer for- 
mulation for the conductance and introducing the DP 
model for decoherent transport. In Section III we de- 
scribe the application of this model to the BL structure 
of PAni and we show in Section IV the effective trans- 
mission results and associated currents in the non-linear 
regime. In section V we present the main conclusion: at 
room temperature decoherence can account for high con- 
ductance in the strongly disordered bipolaronic lattice. 
The mechanisms that can contribute to decoherence are 
discussed in an Appendix. There, we asses the possible 
role of interchain coupling and various forms of electron- 
phonon interactions. Moreover, a detailed treatment of 
these specific models allows us to enlight their differences 
and similarities with the well-known Marcus-Hush theory 
1431 for vibration-assisted electron transfer. 



II. EFFECT OF DECOHERENCE IN THE 
CONDUCTANCE 

The Landauer formulation requires to calculate the 
transmission of the system. Assuming that the sample's 
Hamiltonian is known one must incorporate explicitly 
the leads connecting to the electrodes. Besides, to take 
account of decoherent interactions with an environment 
DP included several phase-breaking fictitious probes [i^ . 
Thus, an original molecular orbital Hamiltonian with N 
orbitals, 



Hence, P^ is an energy uncertainty associated to a de- 
cay rate of the local state at site / described by the PGR. 
We drop any possible dependence on I simplifying the de- 
scription. Since / is the identity operator, 77 can be taken 
as an infinitesimal imaginary part of the local energy, 
Ei — !■ Ei — IT], resulting in a decay to the environment in 
the same sense as the P's above. 

Given the effective Hamiltonian we have the retarded 
and advanced Green's functions in terms of the real en- 
ergy variable s: 
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(6) 
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where Hoff . may be non-linear in e and non-hermitian for 
Tj ^ 0. The imaginary part of the Green's function enable 
to evaluate the local density of states at site i as: 



N,{e) = -- lini ImGfi(£). 

TT ))->-0+ 



(8) 



In a closed system, this is also obtained diagonalizing the 
Hamiltonian. 

According to the optical theorem, the local density of 
states, 
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becomes an effective Hamiltonian that incorporates the 
leads and the interactions with the environment: 

i?cff. = [Hs - ir; /) + El + + E^, (2) 

where El = Elc^Cj^ and E/j = Urc^cn are, respec- 
tively, the self-energies operators describing the escape 
to the left and the right leads obtained through a Dyson 
equation. 



{Eo - IT]) - EL(fl) 

(e) -irL(m(e), 



L{R) 



(3) 
(4) 



where ^l{r) results proportional to the escape rate, and 
hence to the Fermi velocity, at the L{R) lead. We include 
the effects of the incoherent processes in the Hamiltonian, 
e.g. electron-phonon or through space electron-electron, 
simply through an imaginary correction to selected site 
eigenenergies 



E 



iPrf,C,+ Q. 



(5) 



can be cast as a fiux between the asymptotic states de- 
scribed by a generalized Fisher-Lee formula. We adopt 
the notation Pl = P^i, P_r = Paat, P^ = P^i, to empha- 
size that each site can decay through different processes, 
e.g. a, /3 e {L, R, cj)} are the possible decay processes tak- 
ing place at sites i,j e {1, N}. With this notation the 
generalized transmission probability results: 

Tc,^M^) = 2P/3,(£)G« (e)2P„,(£)G5(e). (10) 

Within the DP model current conservation is imposed 
at each site and energy e. This requires an incoherent 
re-injection of every electron that has decayed through 
the P^j. This leads to the evaluation of the kernel W~^, 
describing incoherent density propagation, where (isj : 

N 

^hj = [(E '^'l>k,'l>t)^i,3 + - hj)]- (11) 

k=l 

From this we obtain the effective transmission through 
the sample: 

N 

Trl = Trn,li + E [W^"^] T^j,Li- (12) 

The right hand side of Eq. [T^ contains two contribu- 
tions: the first one represents electrons that propagate 
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FIG. 2. Diagrammatic representation of benzenoid rings dec- 
imated to obtain equivalent renormalized units in one dimen- 
sion. Incoherent channels of D'Amato-Pastawski model also 
shown. 

coherently through the sample, the second term contains 
the incoherent contributions due to electrons that suffer 
their first collision at site i and their last collision at site 

i- 

Since, T(e) piles up all vertical processes into the 
energy e, one can finally approximate the net current 
through the sample by: 

l^^ ffnL{e)[fL{e)-Me)]de (13) 



—TRL{eF)V = GV, 



(14) 



where fii(L) is the Fermi distribution function at the right 
(left) lead and the factor 2 accounts for spin degeneracy. 
Here, we assumed that an electron with energy e coming 
from the left lead arrives to the right lead with the same 
energy, which is only true for very small F^. The second 
line is the linear approximation for infinitesimal voltage 
V and temperature, where G accounts for conductance. 



III. THE PANI MODEL 

We consider a fully protonated BF, which we expect 
to correspond to the highly conducting emeraldine salt. 
By decimation of the benzenoid rings as it is shown in 
figureHl we reduce the PAni emeraldine salt chain to one- 
dimensional effective system [3l|. Each ring is replaced 
by the proper renormalized sites at the place of the para- 
Carbon atoms. 

The sample Hamiltonian results: 



When j = 3s and j = 3s — 1 we have the renormalized 
parameters for para-C p^-orbitals: 



VooVl 



{e-Eo){e-Eo - 



En 



and 



Ei — E.„ 



s-Eo- 



(16) 



(17) 



Eo 



where Eg and Ep are bare site energies for electrons in 
the Pz-orbitals of ortho-C and para-C respectively; Voo is 
the hopping between ortho-C and Vpo is the hopping be- 
tween a para-C and ortho-C. In this work we use the tight 
binding parametrization of Vignolo et al. (isj for the 
bipolaron lattice model of base-emeraldine doped with 
HCl. 

We will consider decoherent sources on effective elec- 
tron p^-orbitals sites by including a constant imaginary 
correction to the site energy as in Eq. [Sj This is the 
most convenient choice for computational purposes. A 
first principles calculation of this imaginary correction 
is beyond the scope of this work. Its complexity lying 
on the multiple effects that must be considered. In the 
appendix we discuss in certain detail two mechanisms: 
interchain tunneling and the effect of torsional modes on 
the crucial tt bonds. However, it is enough to resort to 
dimensional arguments based on spontaneous symmetry 
breaking of the quantum coherent state (4^ . The reason- 
ing takes into account that quantum bosonic modes with 
energies < ksT should be occupied by many quanta, in- 



dicating that r0 should be of the order of ksT [2j|. In 
accordance with this general framework, Rebentrost et 
al. obtain a comparable estimation for F^ for excitons 
in photosynthetic complexes interacting with a phonon 
bath [4^. Indeed, experiments on DNA strands fits F^ 
of this order using the DP model [i^. Thus, we fix an 
effective F^ on each effective site such that the energy 
uncertainty per orbital site is F^ — ksT. We will see in 
the next section that small variations of the precise value 
of F^ have little impact on conductance. 

Right and left leads are described by Eq. |3] choosing 
Eq — and V — 5eV to observe the appropriate band- 
width of interest [— lOeV, lOeV]. Furthermore, we calcu- 
lated the Fermi energy by diagonalizing the exact tight- 
binding Hamiltonian for different configurations. For all 
possible chain arrangements, the Fermi level is nearly the 
same, around the average. 



N 

Hs^Y. i^^^t^o + V,c+cj+, + Vjc+^,c,) . (15) 

When j = 3s + 1 with s positive integer, Ej is the Nitro- 
gen p^-orbital energy and Vj is the tt binding energy (hop- 
ping) between the nitrogen and the para-C p^-orbitals. 



IV. NUMERICAL RESULTS 

We have performed a detailed analysis of the conduc- 
tance properties of the BL model of polyaniline emeral- 
dine salt. Due to the fact that, according to experimen- 
tal data, PAni chains seem to have an average length 
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FIG. 3. Conductance for both: (a) T=OK, and (b) T=300K 
according to DP model for a 400 rings long PAni-HCl chain. 
We also show the T=OK Density of state for comparison pur- 
poses. 
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FIG. 4. Fluctuations in total conductance at T=300K in 
the main peak around the Fermi energy with the fraction of 
quinoid rings along the chain. 




of 400 rings [47| , we have taken that number in our nu- 
merical calculations. However, it should be noted that 
our results do not depend critically on this parameter. 
We first calculated the coherent transmission probabil- 
ity as a function of energy for a set of chain configura- 
tions drawn from the representative ensemble. Results 
are in full agreement with the those of Farchioni et al. 
[T3 . [isl lisj and evidence the mobility edges induced by 
correlated disorder in this 1-D system [9J . While there is 
an appreciable density of states at the Fermi energy, it 
corresponds to locahzed states. Indeed, according to Fig. 
[3]-a, the Fermi energy is far away from the extended state 
region. Fig. [3]-b shows the drastic differences in conduc- 
tance once that decoherent processes are taken into ac- 
count. Conductance at the Fermi energy now becomes 
appreciable for any configuration. These results show 
that metallic transport is possible within a purely BL 
model through an environment-assisted transport psl j. 
Even within a model of perfectly conducting PL islands 
bridged by BL strands, the calculated chains can be taken 
as representative of such transport bottleneck. One ex- 
pects that small differences in quinoid ring concentrations 
would appear due to natural fluctuations on the oxida- 
tion degree previous to doping. In Fig. 2] we show the 



FIG. 5. Conductance for a 400 rings long PAni-HCI chain as 
a function of F^. The value for F^ at T = 300 A" also shown 
in dash line. 



resultant conductance for various quinoid concentrations 
and found no signiflcant changes in transport. 

We also studied the behavior of total conductance as 
a function of decoherence rate. In accordance with re- 
cent works [1^ H^, three regimes can be appreciated. 
Starting from T=OK, as the temperature rises dephasing 
events become more successful in the destruction of lo- 
calization caused by coherent interference at the Fermi 
energy, rising the total conductance of the system. In 
this regime, transport rate increases as the energy uncer- 
tainty associated with temperature is increased. How- 
ever, there is a F^ value for which the conductance is 
maximal. If the temperature is increased further, the 
associated energy uncertainty becomes larger than the 
terms of the system Hamiltonian (characteristic hopping 
and site energies), and the decoherent process now are 
able to suppress transport. This is commonly known in 
the literature as the quantum Zeno effect. In Fig. [S]we 
show our results for the dependence of total conductance 



6 




200 400 600 800 1000 

Ring Number 



FIG. 6. Total conductance at T=300K as a function of the 
chain large. As an inset, it is shown the scale difference be- 
tween coherent and total conductance. 



with r^, in which the three regimes above described are 
clearly seen. The thermal energy, in this case, lies on an 
area of great influence on the total conductance of this 
system, and therefore decoherent processes should not be 
neglected. At room temperatures, the PAni BL is safely 
placed in the range of thermally assisted transport, and it 
is clear from the figure that small variations in the exact 
value of do not alter the outcome significantly (note 
the log- log scale). 

We studied the dependence of conductance on the 
chain length. The results are shown in Fig. |6l With 
the exception of some fluctuations at short lengths, con- 
ductance at T=300K decreases as the reciprocal of the 
chain length, as expected for an Ohmic system. The fit- 
ting gives G/ {2e^/h) = 1/ (20, 67Vfl,), where Nr is the 
number of rings of the chain. The log scale in the in- 
set figure emphasizes the drastic difference between the 
full conductance and that restricted to coherent tunnel- 
ing processes. The coherent conductance decays expo- 
nentiall y a s expected for a one dimensional disordered 
system [50[. In our case the localization length is small, 
G/ (2e2//i) = 22 X We'^-^^^'^, which implies that even 
for very short disordered polymer chains, transport does 
not take place unless decoherence processes are involved. 
Indeed, conductance decays a factor 1/3 for every two 
rings. Therefore, in a model of islands as that mentioned 
in the introduction, destruction of localization by deco- 
herence would have a fundamental role. However, our 
results go further and evidence that even a fully BL PAni 
would sustain strong electronic transport. 

Fig. [7] shows the net current through 50 random chain 
configurations of bipolaronic PAni. The current was eval- 
uated by using Eq. [T2] for symmetric voltages at room 
temperature. As can be seen, the behavior of different 
chain configurations is very similar and close to the lin- 
ear regime. All chain configurations exhibit an apprecia- 
ble conductance with an average value of 7, 7xl0~^ri~^. 
Taking the PAni-HCl density as l,329g/cm3 [5l| we get 




Voltage (V) 

FIG. 7. Average current (black line) and currents for 50 dif- 
ferent configurations (shadow lines) are shown. 

a conductivity of cr = Slfi^^cm^^. This result is slightly 
higher than the experimental ones, lf2~^cTO~^ ^ cr ^ 
2QQ.~^cm~^ (S^l which is reasonable because in the cal- 
culation of the conductivity we are considering the con- 
ductance of the ideal case of chains directly connected to 
the leads. 



V. CONCLUDING REMARKS 

In this work we have discussed electron conductance 
in a doped PAni. We show that the PAni ground state 
configuration, the BL, has high conductance even in pres- 
ence of disorder provided that decoherent processes are 
included. This is done without leaving the convenient a 
la Landauer approach by using the generalization intro- 
duced by D'Amato and Pastawski [32i| where an effective 
transmission accounts for decoherent processes. While 
our formulation accepts further improvements, it pro- 
vides an answer from the robust description of Keldysh 
formalism within a minimal parametrization. Roughly 
speaking, decoherent processes split each chain into a se- 
ries of portions whose length is given by the decoherence 
length |53] ■ These define the elemental conductivities 
from which the sample's Ohmic transport builds on. 

For many years, it has been assumed that conduction 
of polyanilines is inseparably linked to the existence of 
a polaronic crystalline structure. However, although our 
main intention is qualitative, we showed that decoherent 
processes are able to give appreciable metallic conduc- 
tion in the more entropically favorable bipolaronic lat- 
tice. For this system, the uncertainty of energy associ- 
ated with thermal processes cannot be neglected in the 
study of conductance, since ksT falls in a region in which 
the interplay between incoherent and coherent dynamics 
results in an increased efficiency of electron transport. 
One might then speculate that only when the thermal 
energy scale becomes smaller than the Coulomb energy 
of the localized states, one would actually start to notice 
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a qualitative difference with an ideal 1-D metal. 

The robustness of the results obtained is evident by 
noting that they neither depend on variations in the ox- 
idation degree of PAni prior to the doping process, nor 
on the particular arrangement of quinoid rings along the 
chain, or on the exact value of the energy uncertainty 
associated to F^. This justifies the fact that good con- 
ducting properties do not depend much on the purity 
of the emeraldine base so that small displacements to- 
ward the leucoemeraldine or pernigraniline are accept- 
able. The evaluation presented in the appendix show 
that, even when interchain coupling can contribute ap- 
preciably to conductivity, the coupling between the pz 
bonds with torsional degrees of freedom is strong enough 
to provide almost all the required decoherence. This hy- 
pothesis seems consistent with the experiments that show 
that adding residues that restrict the torsional motion 
would also diminish the conductivity as compared with 
the unmodified bipolaronic lattice |54| . 

We do not attempt to rule out the presence of phase 
segregation into metallic polaronic islands and "insulat- 
ing" bipolaronic domains. However, these last strands 
constitute the bottle-neck where thermal decoherent pro- 
cesses activate the conductivity. Moreover, our results go 
further ahead and evidence that bipolaronic chains can 
sustain electronic transport by themselves. In fact, based 
in our simulations we can estimate bulk conductivity for 
these chains and arrive to a remarkably good value as 
compared with experimental data. 
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VII. APPENDIX 



A. Interchain hopping 

We start considering the effect of Vx, an interchain 
hopping at site j. Any neighboring chain can act as an 
"environment" for an electron at this site. This is be- 
cause an electron jumping into a side chain (see Fig. |8]) 
has two options: 1) to escape towards this alternative 
propagation channel and never return. This is obviously 
decoherent as it can not interfere any longer with the 
main pathway [s^ [s^. 2) to return after having an er- 
godic walk on the side chain. In this case it is just the 
excessive amount of interferences and anti-resonances in- 
volved that leads to a decoherent description 'ss']. Each 
node in the plot corresponds to a multi-chain electronic 
state. Notice that the interaction structure looks the 
same as in the local phonon picture discussed in Ref. 
[25I [26| . Using Eq. |2l we have for self-energy de- 
scribing this coupling: 



\Vx? 



U; 



X 



(18) 



where Ej, Vjj^iaie site and hopping strength within the 

chain. Thus the interchain rate can be expressed as F^ = 
2 

Tj where Tj and the imaginary part of total self 

energies at site j. We may evaluate an estimate for the 
typical Fj by disregarding localization and considering 
that the side chain is an infinite PAni strand and using 
the representative values of i? ~ ~0.3eV for and site 
energy and V ~ —3.6eV for intrachain tt bonds. Thus, 



(19) 



We might wonder which range of values would be re- 
quired from Vx to yield an energy uncertainty of the 
order F^ ~ ksTji where kBTn stands for room tem- 
perature energy. The use of the discussed values yields 
Vx < V/12. While we can not ensure that this is the 
case, in every site of the PAni chain, uncertainty energy 
associated with interchain coupling is not too far below 
thermal energy and, therefore, it is it is not negligible. In- 
deed, there could be sites in which interchain couplings 
are stronger, and they would contribute substantially to 
the decay of local states. 




In this section we attempt to elucidate the meaning 
of F^ by presenting two physically meaningful sources 
of decoherence which must be present in our system: 
interchain coupling and a simple but general model of 
electron-phonon interaction. These presentation also al- 
lows to compare the similarities and differences between 
the decoherence rate and the electron-transfer rate in the 
Marcus-Hush model. 



B. Torsional Phonon coupling 

Certainly, vibrational degrees of motion are natural 
sources of decoherence. We analyze them by introduc- 
ing a simple model for electron-phonon couplings that 
enables the evaluation of the corresponding contribution 
to F^. From the geometrical inspection of the molec- 
ular structure, it is obvious that torsional strains on 



\J,S 



\j-l,n+2) \j,n+2) \j +l,n+2) 




l,/+U'> 



FIG. 8. Interchain hopping at site j. States are written in 
Dirac notation including quantum numbers s which label dif- 
ferent PAni chains. This representation illustrates the simi- 
larity with Fock-space representation of the electron-phonon 
system. 



benzenoid rings disrupt tt bonds between orbitals 
of para-Carbons and Nitrogens. Their overlap depends 
on the angle 9 between the orbital axes. As a re- 
sult, the corrected hopping energies can be written as 
V = V°cos{e) ~ V°{1 - 612/2). The natural frequency ujg 
of this torsional motion determine the vibrational energy 
of benzenoid rings. A self consistent description requires 
that the restoring force lujgO, written in terms of the mo- 
ment of inertia / of the benzenoid ring, should coincide 
with the net change in the electronic energy described by 
the tight-binding model. In this case it yields 

leading to hug ~ 2 x lO^^eV < fc^TR. 

In terms of the second quantization operators b — 
y/luJ^j2h{e + ie/oje) and b+ = ^Iuj0/2h{e - iO/ojee) we 
get the perturbation given by the coupling Hamiltonian: 



'el-ph = - jhwg (b+ + i)^ (Sj'^-i + Sj'j-i) 



(20) 



A Fock-space representation of this interaction Hamil- 
tonian is represented in Fig. [S] Notice the similari- 
ties and differences with the representation of the lin- 
ear electron-phonon interaction discussed in Ref . [1^ [2^ 
and the interchain coupling. In the present case, the ef- 
fect of the perturbation on the state on a local site j can 
evaluated with the FGR: 



2^ , I „ 2 



5[{£ + nhLUe) - (Ej, +n'huje)] (21) 

where \j,n) — ( ) 1'^'* ^^"^^ \^) ^^e electron 
and phonon vacuum and n label the number of vibra- 
tional quantums whose thermal probability is P{n). In 




i-\,n-2) \j,n-2) \j +U-2) 



FIG. 9. Fock-space representation of state | j, n) and its sur- 
roundings. The middle row represents electronic states with n 
phonons in the PAni chain. Lower and upper rows represent 
the same chain but with different numbers of phonons. Black 
dotted lines are electron-phonon couplings. 



the case of interest, we consider electrons at the Fermi 
level, Ep- Thus, after energy integration and using the 
thermal average {{n)) = n = '^P{n)n for the expecta- 
tion number of n, the decay rate becomes: 



- = Tk(f^o?{{n^ + 4n + 2)N{Ef - 2hujg) 
Tj Ion 

+ 2n^N{EF + 2hujg) + {8n^ + 8n + 1)N{Ef)} (22) 

We must highlight that the quadratic dependence on 
displacement in the electron-phonon interaction is the 
responsible for the selection rules that can be appreciated 
in the last equation. Electrons are allowed to interact 
with environment only by absorbing or emitting phonon 
pairs. This is shown in Fig. [9l However, without much 
loss of generality that fc^T ^ huig, so it is possible to 
approximate Ep ~ Ep ± 2hujg and n « jS^- As a result, 



- = ^{HuJerNiEp) 
Tj 8 ft 



12 



12 



(23) 

The evaluation of the corresponding F^ becomes triv- 
ial in this high-temperature regime. 



h 1 Stt , , , 

-- - —N{EF){kBT)\ 



(24) 



Here, it is crucial to notice that for highly localized 
states the imaginary self energy results mainly from the 
decoherent process described above. Thus 



1 



N{e) « - 



^ (e - EoY 



(25) 
(26) 
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FIG. 10. Fock-space representation of Eq. [30] and its corresponding semiclassical picture in terms of vibronic degree of freedom 
X as used in the Marcus-Hush model for electron transfer reactions. The potential surface at the right becomes shifted when 
represented in the polaronic basis allowing to define the reorganization energy Er- In this model, the electronic coupling Vab 
must be small enough to remain perturbative. Left panel is also used to represent the decoherence by electron-phonon coupling 
when Vg < Vab- 



Therefore, from Eq. [24l 



(27) 



Thus, any low frequency modes yielding a quadratic 
dependence of the electronic energy on the displacement, 
which can be more general than expected, leads to an im- 
portant consequence: for localized regime, it would tend 
to provide an energy uncertainty (decoherence) of the 
order of the thermal energy. A similar behavior remains 
valid if one relaxes the localization requirement to that of 
a sharply peaked resonance. A simple example is a sharp 
resonance in a one dimensional system js?} in which the 
local density of states could be written: 



N(e) 



-Im 



1 



1 



n e — Ej + ir^ — q;(A — iP, 



bulk) 



(28) 



where Tbuik represent the escape to the rest of the tight- 
binding chain and A gives the energy shift due the pres- 
ence of the other sites. For e within the band edges, this 
equation results: 



N{e) 



1 



1 + 2arb„(fc/r 



ttF 



( ^"%"'"^ ) + (1 + 2ar,„,fe/r)^ 



(29) 



which in the limit of large F compared with aTi,uik, gives 
N(e) « I/ttF. This limit is achieved at room temperature 
whenever a ^ 1. 



C. Comparison to Marcus-Hush theory 

In the previous section we deemed with the physi- 
cally relevant situation of a quadratic interaction with 



the vibrational coordinate. This is not conceptually dif- 
ferent with the standard linear electron-phonon coupling 
used to describe the Franck-Condon effect [H^ and the 
electron-transfer process [5i|. All these physical pro- 
cesses are contained in a simple Hamiltonian 



H = E Ejc+Cj +hLOo{b+b+\] -Vg[b+ + b 

J=A,B V / V 



VAB{c\cg + c]^l 



B'^B 



(30) 



whose interactions in the Fock space are represented in 
Fig. [TUl The electron transfer problem is best repre- 
sented resorting to the polaronic transformation which 
would diagonalize the Hamiltonian but for the tunnel- 
ing described by Vab- The essence of an estimation 
of the electron transfer rate is a FGR evaluation of 
the tunneling between the electronic states A and B in 
the regime of weak coupling non-adiabatic limit fnjo <C 

kBT,\VAB\ < \Vg\- 



'^-^\Vab\^F{/^E)] (31) 



TA- 



where AE = Ea - Eb, Eb = Eb - V^/Hujq and F{AE) 
is a density of directly connected states denominated 
Franck-Condon factor. Thus it satisfies 



F{AE)dAE ^ 1. 



(32) 



F{AE) is estimated resorting to a thermal average and 
following the Marcus original treatment which interprets 
the transition probability according to a Landau-Zener 
formula. Thus 



F{AE) 



1 



y/AT:ErkBT 



exp 



{AE - Er) 
AErkBT 



(33) 



10 



where the reorganization energy Er is indicated in the 
plot. 

In contrast to this treatment, in a decoherence problem 
one focus on estimating how the electron-phonon interac- 
tion degrades the standard coherent Rabi oscillation [60j. 
This describes an electron jumping forth and back be- 
tween states B and A and atenuates within a decoher- 
ence time T0. Similarly to Eq. (3.14) of Ref. t33|] a FGR 
evaluation gives: 



^='-i\Vf{{{n+l)N{EB 
~<t> 



huJo)+ 7iN{Ee 



■N{Eb] 



hwo))) 

(34) 



where (()) stands for thermal average. The approxima- 
tion involves a high temperature limit and again, the 
square brackets indicate a density of directly connected 
states. As in previous section, the assumption that 



phonon induced electron energy uncertainty leads to the 
self-consistent condition of Eq. 



h 

2^ 



Wa 



2TTkBT 



(35) 



which is valid provided that hajQ <C fc^T, |T/g| <C Vab- 
Notice that in an ab-intio parametrization of the tight- 
binding Hamiltonian, the coupling constant results from 
evaluating the dependence of the parameter on the ap- 
propriate generalized coordinate, e.g. Vg = dEs /d6. 

Thus, in both problems, electron transfer in presence 
of a some reorganization energy and electron transport 
with decoherence from a phonon bath, the Hamiltonian 
is the same. However, since the calculated observables 
are different, the term used as perturbation in the FGR 
differ. In the first case the perturbation is the electron 
jump Vab, while in the decoherent situation the pertur- 
bation is the electron-phonon coupling constant Vg. 
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